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Abstract 

An upgraded analysis of e, Xd and e'/e, using the latest determi- 
nations of the relevant experimental and theoretical parameters, is 
presented. Using the recent determination of the top quark mass, 
mt = (174 lb 17) GeV, our best estimate is e'/e = 3.1 ± 2.5, which lies 
in the range given by E731. We describe our determination of e'/e and 
make a comparison with other similar studies. A detailed discussion 
of the matching of the full theory to the effective Hamiltonian, written 
in terms of lattice operators, is also given. 
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1 Introduction 



The understanding of mixing and CP violation in hadronic systems consti- 
tutes one of the crucial tests of the Standard Model. In the last few years 
considerable theoretical and experimental effort has been invested in this 
sector. 

On the experimental side, more accurate measurements of the mixing 
angles are now available and the mass of the top quark, recently discovered, 
is constrained within tight limits Still, in spite of very accurate measure- 
ments, the experimental results for the CP violating parameter e'/e are far 
from conclusive [0, 

On the theoretical side, the complete next-to-leading expressions of the 
relevant effective AS" = 1, AS" = 2, AB = 1 and AB = 2 Hamiltonians have 
been computed 0^0; thus reducing the theoretical uncertainties^ More- 
over, there is now increasing theoretical evidence that the value of the pseu- 
doscalar 5-meson decay constant is large, J'b ~ 200 MeV, and that the 
B^-B^ mixing parameter Bb is quite close to one. Still, the evaluation of 
other hadronic matrix elements is subject to large uncertainties, which are 
particularly severe for e'/e, where important cancellations of different contri- 
butions occur, for large values of the top mass. A real improvement in the 
calculation of the hadronic matrix elements, from lattice simulations, or with 
other non-perturbative techniques, is still missing. In spite of this, we believe 
that, given the experimental and theoretical novelties, it is time to present an 
upgraded analysis of mixing and CP violation for K- and 5-mesons, along 



the lines followed in refs. [|1^, We will also have the opportunity of pre- 
senting our results for the Wilson coefficients of the operators appearing in 
the AS = 1 Hamiltonian in different regularization/renormalization schemes, 
which may be useful for further phenomenological analyses by other authors. 
We also intend to present and clarify several issues related to the matching 
conditions of effective Hamiltonians to the full theory, to the use of the lattice 
regularization, to the value of the mass of the strange quark and to quark 
thresholds in the evolution of the Wilson coefficients. Particular emphasis is 
devoted to a realistic evaluation of the uncertainties in our predictions. 



^ Indeed only the top contribution to the AS" — 2 Hamiltonian is fully known, at the 
next-to-leading order. To our knowledge, the other contributions have been only partially 
computed fel. 
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The results given here contain several improvements with respect to our 
previous work on the same subject The main differences are the follow- 
ing: 



Upgraded values of the experimental parameters entering in the phe- 
nomeno logical analysis, such as the 5- meson lifetime tb, the B^-B^ 
mixing parameter Xd, the CKM matrix elements, {\Vcb\, \Vub\/\Vcb\), 
etc., have been used. 

The value of the strange quark mass has been taken from lattice 
calculations, thus making a more consistent use of lattice results for 
the 5-parameters of the relevant penguin operators |P^|. 



The relevant formulae, necessary to express the effective Hamiltonian, 
derived in the full theory, are given in terms of lattice operators. 

Assumptions and relations which have been used to evaluate the B- 
parameters of the operators entering the calculation, in particular those 
that have not been computed either on the lattice, or with other non- 
perturbative techniques, are critically reviewed. 

Uncertainties, coming from the matching between lattice and conti- 
nuum operators, are discussed, together with those due to the choice 
of the renormalization scale, Aqcd, the mass of the top quark, etc. We 
also discuss different determinations of e'/e originating from the choice 
of different continuum regularization schemes. 

All the results are presented with an estimate of the corresponding 
errors. These errors come from the limited precision of the measured 
quantities, e.g. tb, and from theoretical uncertainties, e.g. the value of 
hadronic matrix elements. 



The phenomenological results of this study have been partially reported else- 
where [p!3| , In this paper we put particular emphasis on some theoretical 
aspects and subtleties related to the construction of the effective Hamil- 
tonian. The reader will also find a practical tool to combine the Wilson 
coefficients with his own preferred results for the operator matrix elements, 
sections 
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The plan of the paper is the following. We summarize in section |^ the 
main results obtained from a combined analysis of e, Xd and e'/e together 
with our estimate of the errors. All the details of the analysis will be given 
after this section. In section ^, we give the relevant formulae of the effective 
Hamiltonians, which control e, B^-B^ transitions and e'/e, and define the 
relevant operators. In section |, a summary of the calculation of the Wil- 
son coefficients of the operators, at the next-to-leading order, is given. The 
matching of the continuum theory with the lattice is addressed in section 
^. The formulae presented in this section have been used to derive the re- 
lation between the mass of the quarks on the lattice and in the continuum. 
In section ^ we give details on the Regularization Independent {RI) renor- 
malization prescription, together with the numerical values of the Wilson 
coefficients in different renormalization schemes. In section |^, the values and 
errors of the 5-parameters, used in this analysis, are given and compared 
with those used in ref. A discussion of the assumptions and relations 

among the i?-parameters is also presented in the same section. The main 
uncertainties of the phenomenological analysis are reviewed in sec. |^. 



2 Results 

In this section, the main results of the present study are summarized. These 
results have been obtained by varying the experimental quantities, e.g. the 
value of the top mass rrit, tb, etc., and the theoretical parameters, e.g. the 
kaon i?-parameter Bk or the strange quark mass mf^^{fj,), according to their 
errors. Values and errors of the input quantities, used in the following, are re- 
ported in tables ^ For the experimental quantities, we assume a Gaussian 
distribution and, for the theoretical ones, a flat distribution (with a width 
of 2 a). The only exception is m^^{^), taken from quenched lattice QCD 
calculations, for which we have assumed a Gaussian distribution, according 
to the results of ref. fl^ . 

Each theoretical prediction (cos 5, sin 2/3, e'/e, etc.), depends on several 
input parameters, which fluctuate independently. In many cases, this pro- 
duces a pseudo-Gaussian distribution of values, see for example e'/e in flg. 
m. From the width of the pseudo-Gaussian, we estimate the errors of our 
predictions. Non-Gaussian distributions may also occur, for example when a 
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Parameters 


Values 


m,(2 GeV) 

^^QCD 

Kfc = AX' 
\Vub/Vcb\ = Xa 
Tb 

{fBBy^)th 

VLib 


(174 ± 17) GeV 
(128 ± 18) MeV 
(230 ± 80) MeV 
0.040 ± 0.006 
0.080 ±0.015 
(1.49 ±0.12) X 10-12 
0.685 ±0.076 
(200 ± 40) MeV 
0.25 ±0.10 



Table 1: Values and errors of the parameters used in the numerical analysis. 



single input parameter dominates the result. In these cases, we still estimate 
the error from the variance of the distribution. 

Using the numbers given in the tables and the formulae given in the 
forthcoming section, we have obtained the following results: 

a) The comparison of the experimental value of e with its theoretical 
prediction will in general correspond to two possible solutions for the value of 
cos 6, see for example ref . |jTO| . Because of the errors of the input parameters. 



the two solutions appear as a distribution of values, with a twin-peak shape. 
The distribution for cos 5 is given in fig. |l|. As already noticed in refs. ||10|, |ll| 
and [|TB], |T^, a large value of /s, combined with a large m^, favours cos 5 > 0. 



When the condition 160 MeV < /sS^^^ < 240 MeV is imposed (/^-cut), 
most of the negative solutions disappear, giving the dashed histogram of fig. 
|I|, from which we estimate 

cos 5 = 0.47 ±0.32 . (1) 

b) cos 5 is correlated to the values of the Wolfenstein parameters p and 
rj. In fig. 0, the contour plot in the p-rj plane is given, with and without the 

/fi-CUt. 
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Constants 


Values 


(jr p 


i.iooo4 X iU UeV 


rric 


1.5 GeV 


rUb 


4.5 GeV 




Sn 2 OeV 




140 MeV 


Mk 


490 MeV 


Mb 


5.278 GeV 


AMk 


3.521 X 10-12 MeV 




132 MeV 


h 


160 MeV 


A = sin 6 c 


0.221 


^exp 


2.268 X 10-3 


ReAo 


2.7 X 10-^ GeV 


UJ 


0.045 




2 GeV 



Table 2: Constants used in the numerical analysis. 



Bk 


^^3/2) 


^1-2 


-83,4 


-85,6 


r(1/2) 


r(3/2) 


0.75 ±0.15 


0.62 ±0.10 


- 0.15W 


1 - 6^*) 


1.0 ±0.2 


1(*) 


1.0 ±0.2 



Table 3: Values of the B -parameters, for operators renormalized at the scale 
H = 2 GeV. The only exception is Bk, which is the renormalization group- 
invariant B -parameter, defined in ref. jl^/; Bg^'^{fi) has been taken equal to 
Bk{i^), at any renormalization scale jj,. In the table, we give Bg^'^{n = 2 
GeV). Entries with a are educated guesses; the others are taken from 
lattice QCD calculations /|7^-/p^/. 5*66 section^ for more details. 
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m,= (174±17) GeV A'**ocd = (330 ± 1 00) MeV 
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Figure 1: Distributions of values for cos 6, sin 2/5, e'/e and e' / e{A^ri)~^ , 
for rrit = (174 ± 17) GeV, and using the values of the parameters given in 
tables 1^. The solid histograms are obtained without using the information 
coming from B^-Bfi mixing, the dashed ones after imposing that 160 MeV 
< fsB]^"^ < 240 MeV. 
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m,= (174±17) GeV A'**ocd = (330 ± 1 00) MeV 
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Figure 2: Contour plots in the plane. The solid, dashed and dotted 
contours contain 5%, 68% and 95% of the generated events respectively. The 
contours are given by excluding or including the fs-cut. Similar results can 
he found in refs. j^, [7^/ . 
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c) The value of sin 2/3 depends on cos 6 

2o"sin5(l — crcos(5) 

sin 2/5 = , ; _ /. 2 

1 + cr^ — 2a coso 

The distribution of sin 2/5 is shown in fig. |T|, without (sohd) and with 
(dashed) the /^-cut. When the /s-cut is imposed, one gets larger values of 
sin 2/? [|TU|. From the dashed distribution, we obtain 

sin 2/5 = 0.65 ±0.12 . (3) 

d) In fig. H, several pieces of information on e'/ e are provided. Lego-plots 
of the distribution of the generated events in the e'/ e-cos 6 plane are shown, 
without and with the /^-cut. At the same time, the corresponding contour 
plots are displayed. One notices a very mild dependence of e'/e on cos 5. As 
a consequence, one obtains approximately the same prediction in the two 
cases (see also fig. ||). In the HV scheme ||33| the results are 

e'/e = (2.3 ± 2.1) x 10"^ No-cut, (4) 

and 

e'/e = (2.8 ± 2.4) x 10"^ /s-cut, (5) 
whereas in the NDR scheme they become 

e'/e = (2.8 ± 2.2) x 10"^ No-cut, (6) 

and 

e'/e = (3.4 ± 2.5) x 10"^ /s-cut. (7) 

Our best estimate, reported in the abstract, has been obtained by averaging 
the results given in eqs.(^) and (0) and by adding a systematic error, due to 
the difference of the central values in the two schemes, to the final result 

e'/e = (3.1 ± 2.5 ± 0.3) x 10"^ /s-cut. (8) 

e) In view of the rapid evolution of the experimental determination of 
the value of the top mass, expected in the near future, we also give e'/e as a 
function of rrit, in fig. ^. 



By comparing the present results to the analysis of ref. [|rT|, one may 
notice the following: 
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m,= (174±17) GeV A'*'ocd = (330± 1 00) MeV 
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Figure 3: Distributions of the events in the plane e' /e-cos5 without and with 
the fB-cut. The corresponding contour plots are displayed below the lego- 
plots. 
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(GeV) 

Figure 4: e'/e as a function of rrit, obtained by applying the fs-cut. The 
zones delimited by the solid and dashed curves represent the allowed regions 
containing respectively 68% and 95% of the generated events. The region 
between the two double lines is the experimental result coming from E731 
e'/e = (7.4 ±5.9) x 10"^ The range from NA31 /§/, e'/e = (23.0 ±6.5) x 10^^ 
is not given in the figure. 



1. The latest values of iKife/Kfel and \ Vcb\ are sensibly lower than the pre- 
vious ones. Consequently, the separation of the positive-negative cos S 
solutions takes place at larger values of 

2. e'/e ~ iKib/Kbl X iKbP- The actual results for e'/e, however, are not 



very different from those given in ref . |Tl| , because we are now using a 
smaller value of the strange quark mass taken from lattice QCD 
and not from QCD sum rules p4l-p9[, as we did before. One has 



new reasons to doubt of those QCD sum rule determinations, which 



have been found to be affected by large higher order |40] or instanton 
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effects |^T]|. It is reassuring that tlie recent QCD sum rule calculation 
of ref. [4^, apparently free of uncontrolled non-perturbative effects, 



agrees with the lattice calculation of ref. Wl 



3. In spite of several differences, the bulk of our results overlaps with those 
of ref. |]T5|. It is encouraging that theoretical predictions, obtained by 
using different approaches to evaluate the operator matrix elements, 
are in good agreement. 

4. On the basis of the latest analyses, it seems very difficult that e'/e is 
larger than 10 x 10"'^. Theoretically, this may happen by taking the 
matrix elements of the dominant operators, Qq and Qs (see below) very 
different from usually assumed. One possibility, discussed in ref. [p!5|] , 
is to take the corresponding i?-parameters to be i?6 ~ 2 and -Bg ~ 1, 
instead of the values i^g ~ -Bg ~ adopted here and in refs. [1^, |Tl], |1^ . 
To our knowledge, no coherent theoretical approach has found Bq so 
large, and we see no reason to use Bq = 2. 



3 — transitions and e'/e 

In this section, we list the expressions of the effective Hamiltonians, respon- 
sible for AS" = 2, AB = 2 and AS* = 1 transitions, given in terms of the 
relevant operators and their Wilson coefficients. From these Hamiltonians, 
one derives the expressions for e, Xd and e'/e, written as combinations of 
the Wilson coefficients times the 5-parameters of the different operators, i.e. 
their matrix elements. The formulae of this section have been used in our 
study. 

1) The effective Hamiltonian governing the AS = 2 amplitude is given by 
n[f/=' = ^M^id^^sf {xlF{x,) + \lF{xt) + 2AAF(x„ x*)} , 

(9) 

where Gp is the Fermi coupling constant and 7^ = 7'^(1— 75) ; A^'s are related 
to the CKM matrix elements by \q = VqiV*f, where 'i' and '/' are the labels of 
the initial and final states respectively (in the present case i = s and f = d); 
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Xq = m?^/ M^r and the functions F{xi) and F{xi, Xj) are the so-called Inami- 
Lim functions [Q, obtained from the calculation of the basic box-diagram 
and including QCD corrections [H3|; F{xt) is known at the next-to-leading 
order, which has been included in our calculation. From equation (^) we can 
derive the CP violation parameter e 



|e|5=o = a^xA^AV sin 6 Xt) + F{xt) [A^X^l - a cos 5)] - F{xc) } 



(10) 



where 



a 



guimkM^ 



(11) 



Here, AM^ is the mass difference between the two neutral kaon mass eigen- 
states. In eq. (p!0|), p = a cos 6, r] = a sin 5 and A, A, p and rj are the 
parameters of the CKM matrix in the Wolfenstein parametrization 
is the renormalization group-invariant 5-factor, to be discussed in sec. ^ 
Bk takes into account all the possible deviations from the vacuum inser- 



tion approximation in the evaluation of the matrix element (K^Kd'jj^s) 
{Bk = 1 corresponding to the vacuum insertion approximation). 

2) The AB = 2 effective Hamiltonian is given by 



n 



|AB|=2 
eff 



p2 



from which one finds 
AM 



167r^ 



Cb^BbA'xHi + a^- 2acos6)F{xt), 

IVl B 



M^Xiid^^bYFix,) 



(12) 



C 



B 



(13) 



where Bb is the renormalization-invariant 5-parameter, relevant for B-B 
mixing. 

3) The AS* = 1 effective Hamiltonian above the charm threshold is given 

by 



n/AS=l 
^eff 



+ 



rQipf 



1 



Cip) 



CM {Q,{p) - Qlip)) + C2{p) (Q2(/^) - Q^2(/^ 

(14) 
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where Q{fi) = {Qi{fx),Q2ifx), . . .), C{fx) = Csl/i), . . .), A„ = KdK*s 

and similarly Ac and Xf, t = — At/A„ and /i is the renormalization scale of 
the operators Qi. A convenient basis of operators |PB|-pU|, when QCD and 
QED corrections are taken into account, is 



Qi 
Q2 

Q3,5 

Q4,e 

Q7,9 



[Sada){V-A) {upUp) (y-A) 
{Sad/3) (V-A) {Uf3Ua){V~A) 
[Sada){V-A) JK^/S^'/s) (V^A) 

q 

{Sadi3)(V-A) '^{Ql3qa){VTA) 

q 

3 

-{Sada){V-A) '^qiQpQls) {V±A) 



Q 



8,10 



-{s^d 



a"/3j(y-A) 



a)(V±A) 



Ql = iSada)(V-A)iCi3Ci3)(^V-A) 

Q2 = {^adf3)(V~A){Ci3Ca){V~A) (15) 

where the subscript (V ± A) indicates the chiral structure and a and /? are 
colour indices. The sum over the quarks q runs over the active flavours at 
the scale /i. 

From "^^i//"^' we can derive the expression for e' 

gi7r/4 



UJ 



uJ-\lmA2y - (1 - Qib) ImAc 



(16) 



V2 ReAo 

where (Im742)' and ImAo are given by 

-- -^Im {v*Vt,) |- (c,Be + ^C,B,^ Z + (c,B, + ^Ca^s) X+ 



Im An 



,1/2 



/2Y 



X 



V 3 ■ 6 ■ 2 
1/2X (CiBl 



CrBr ( — + - + -)+ CsB;,'^ {2Y + - + -]- 



8^8 



1/2 



X 



6 



CgBg — + 



C2B2 X 



(17) 



and 



(Im A2)' = -G^Im {V*V,,) \c,bT (| - f ) + 
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CsBriY-^)+C,Br'-§]. (18) 

uj = lie A2/Iie Aq = 0.045 and we have introduced {lmA2)' defined as 

lmA2 = {lmA2y + niB{ujlmAo). (19) 
^iB = 0.25 ±0.10 represents the isospin breaking contribution, see for exam- 



ple ref. |£0[- With the subtle cancellations, occurring for a large top mass, 
nit ~ 175 GeV, a precise knowledge of the value of Qib becomes impor- 
tant. For example, by taking Qib = 0, one gets e'/e ~ 5 x 10~^, instead of 
2.8 X 10~^. In view of this, further efforts should be devoted to improve the 
accuracy of the theoretical predictions for isospin-breaking effects. 

The numerical evaluation of e'/e requires the knowledge of the Wilson 
coefficients of the operators and of the corresponding matrix elements. The 
Wilson coefficients have been evaluated at the next-to-leading order, using 
the anomalous dimension matrices given in refs. H]-[pi, and the initial condi- 



tions computed in refs. [0, ^ (and given for HV and NDR in refs. 0-01). 
The matrix elements of the operators have been written in terms of the three 
quantities (see eqs. (O) and ([T8|)) 




K - Ml) , (20) 
Ml V _ /0.15GeV' ^ 



12X — — ^ ' (21) 



l] Y, (22) 



fK 



fn 



and a set {Bi} of i?-parameters (in our normalization /^r = 132 MeV). The 
numerical value of the 5-parameters have been taken from lattice calcula- 
tions [|l^-[^ and multiplied by suitable renormalization factors to take into 
account the difference between HV (NDR) and the lattice regularization 
scheme. For those i?-factors that have not been computed yet on the lattice 
we have used an educated guess, which will be discussed in detail in the 
following. We observe that, in eqs. (|l^ and (0), only nine coefficients {B- 
parameters) appear since we have used the relation Qio = —Qs + Qi + Qg, 
which is valid below the bottom threshold. 
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4 Wilson expansion and renormalization 



This section is devoted to the Wilson coefficients of the operators appearing 
in the effective AS = 1 Hamiltonian, which we have computed at the next- 
to-leading order, including QCD and QED corrections. We do not discuss 
the determination of the coefficients of the operators relevant for e, since this 
point was explained in great detail in ref. |I0|] . 

The effective Hamiltonian is defined by the Wilson Operator Product 
Expansion (OPE) of the T-product of two weak currents: 

{F\n[f='\l) = gl,/8jd'xD^^{x\M'^){F\T{j,{x),Jl{0))\I) 

^ J2a{Mw){F\Qi{Mw)\I). (23) 

i 

— * 

The Wilson coefficients C{Mw) = {Ci{Mw), C2{Mw)---) are found by match- 
ing, at 0{ae) and 0{as), the current-current and penguin diagrams com- 
puted with the W and top propagators to those computed with the local 
four-fermion operators in the effective theory. The C(/i) are expressed in 
terms of C{Mw) through the renormalization evolution matrix M^yj^ 



where 



with 



and 



C{fi) = W[fi,Mw]C{Mw), (24) 



Mw] = M[/i]f/[/i, Mw]M'[Mw], (25) 



MMHi.£^^)(i.t^i)(i.^^) pa) 



At the next-to- leading accuracy, W[fi,M\Y] is regularization scheme de- 
pendent. 

^ We have properly taken into account the beauty threshold in the evolution matrix 
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We now give the basic information necessary to compute the matrices 
defined in eqs. (|25|), (|26|) and (pT]). 

The matrix ?7[/i, Mw] in eq. is given by 



a. (Ml 



(28) 



The matrices P, J and K are solutions of the equations 



J- 



P 



J, 



2/?o 



7, 



{0)T 



7 



(0)T 



^^(0)T _ 7i 



(i)r 



2Pi 



'2 7s 



2Po 



(29) 



(30) 



i^, 7f T = + 7^°)^ J + t^^P + [7i°)^, JP] - 2AP - |p7f ^ . (31) 

The anomalous dimension matrix, which includes gluon and photon correc- 
tions, can be written as 



7 



47r 



7i°) + -T^ir + 



a 



(0) 



tts tte _^ (1) 



(32) 



where each of the 7^V^ is a 10 x 10 matrix. In eqs. (p^)-(|3T|), Pq and Pi 
are the first two coefficients of the /9-function of a^. At the leading order, 
the QCD anomalous dimension matrix, including QCD penguins, has been 



computed in refs. 0. The electroweak anomalous dimension matrix at 
the same order can be found in refs. [Q, |^ and Two groups have 
computed the anomalous dimension matrix at the next-to-leading order, by 
calculating all the current-current and penguin operators at two loops up to 
order a'^t and aeUst 0^@]- AH the details of the computation can be found 
in these references. Here we only give the numerical results for the Wilson 
coefficients in three different regularization/renormalization schemes, which 
are commonly used in the literature, namely HV, NDR and RI: 



1. HV denotes the MS renormalization prescription in dimensional reg- 
ularization, with 75 defined according to the 't Hooft and Veltman 



prescription ||33| 
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2. NDR is the MS scheme in Naive Dimensional Regularization, with 
anticommuting 75; 



3. The RI (Regularization Independent) scheme requires a few words of 
explanation. In its simplest realization, one fixes the renormalization 
conditions of a certain operator, by imposing that suitable Green func- 
tions, computed between external off-shell quark and gluon states, in 
a fixed gauge, coincide with their tree level value. For example, if we 
consider the generic two-quark operator Or = ip'i^ip, we may impose 
the condition 

ZT{p\0r\p)\,2^-,^ = {p\Or\p)o, (33) 

where F is one of the Dirac matrices and {p\Or\p)o is the tree level 
matrix element. The extension to more complicated cases, including 
four-fermion operators and operator mixing, is straightforward and will 
be discussed in detail in sec. |^. The RI procedure defines the same 
renormalized operators, i.e. the same coefficients, in all regulariza- 
tion schemes {HV, NDR, Pauli-Villars, lattice), provided they are 
expressed in terms of the same renormalized strong coupling constant. 
However, the coefficient functions now depend on the external states 
and on the gauge used to impose the renormalization conditions^. Thus 
the external states and the gauge must be specified. To make contact 



with the non-perturbative method proposed in ref. [Q, we will use 
the Landau gauge. The external states will be different for different 
operators, and will be given in sec. ||. 

The RI scheme is particularly convenient for matching the coefficients 
of the operators to the corresponding matrix elements computed in 
lattice simulations, as extensively discussed in ref. The details 

of the implementation of this renormalization scheme will be given in 
sec. where we also present the relations between the renormalized 
operators and coefficients computed with the RI scheme and with the 
usual MS prescriptions. 



■^The "physical" , hadronic matrix elements of the effective Hamiltonian are independent 
of the external states and of the gauge used to fix the renormalization conditions up to 
the order at which one is working, in our case up to the NLO. 
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5 The effective Hamiltonian written in terms 
of lattice operators 



Hadronic matrix elements from lattice QCD have been widely used to predict 
several quantities of phenomenological interest in weak decays, deep inelastic 
scattering and heavy flavour physics. In this section we discuss the renor- 
malization properties of lattice operators and clarify some issues related to 
the construction of physical amplitudes starting from the matrix elements 
computed on the lattice. For definiteness and simplicity, we will refer to 
two-fermion operators and to the four-fermion A J = 3/2 operator Q^, which 
renormalize multiplicatively in the continuum^. The relations derived below 
are, however, valid in the general case. These relations have been used to 



compute the strange and charm quark masses in ref. [|T^ . 

Let us consider first the renormalization of the generic two-quark operator 
ipTip, where F is one of the Dirac matrices. At the next-to-leading order 
(NLO), the generic, forward, two-point Green function, computed between 
quark states of virtuality = /i^, has the formQ 

(0) /,m\2 



IGtt 



167r2 \ 2 V vr y 



2 




where the dots represent terms beyond the next-to-leading order and terms 
of 0(a), which we will assume negligible in the following; Fq is the zeroth- 
order Green function, and g\{a) = 6//3 is the bare lattice coupling con- 
stant; 7^°-* and 7'-^'' are the leading (regularization- independent) and next- 
to-leading (regularization-dependent) order anomalous dimensions, respec- 
tively; A = A(a) is the lattice gauge parameter of the gluon propagator 

^Q'^ is related, via a chiral transformation, to the AS* = 2 operator of the effective 
Hamiltonian, introduced in eq. (^) . We have chosen to discuss to avoid the complica- 
tions present in the calculation of the effective AiS* = 2 Hamiltonian, due to the matching 
conditions in the presence of a heavy top quark. 

^We work on an Euclidean lattice. 
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(^n^i/(g^) = —5^1/ + (1 — \)q^qvl(f^- It obeys the renormalization group 
equation 

' -a'-^--P.i9lia))^-pfM + -^ Pl-'-^^- (35) 



A(a) da 167r2 ' ' 

We have introduced a scale-dependent gauge parameter in order to define 
a gauge-independent anomalous dimension and simplify the renormalization 
group equation for r(/ia), see eq. (|38| ) below. 

The lattice coupling constant obeys the equation 

/alia) _ ^/ 2. XX ^?i(«) , 

where /3o,i are given by 

(lliV-2n.) ^ 34^. 10,, iN"^-!) 

and n/ is the number of fiavours. Equation (34) guarantees that all the matrix 
elements can be made finite (as a — * 0) by multiplying the bare operator by 
a suitable renormalization constant, obtained by fixing the renormalization 
conditions for Op. 

From the above equations, we find 

(a^ - /^(^i(«))^^ - M{a))X-^ - i{gl{a)))T{^,a) = 0, (38) 
with 



-'iAC')) = ^7'°' + ^7'". (39) 



In view of the comparison with some continuum regularization, it is con- 
venient to expand the bare Green function r(yua) of eq. (34) in terms of the 
continuum minimal subtraction {MS) coupling constant, evaluated at the 
scale vr/a. The continuum MS coupling as(7r/a) is related to the lattice bare 
couphng (a) = g1{a)/4:7i by the equation 

^ ' -A + ...), (40) 



(a) as(7r/a) ^ Ati 
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where A is a numerical constant. With this substitution, eq. 



and 



47r 



a. 



(47r)2^ 



7l(«s) = f^7^°^ + 



47r 



2 ' 



(47r 



becomes 
(41) 

(42) 



(43) 



By changing the expansion parameter, one also has to change the two-loop 
anomalous dimension p[-p[, see also eq. (34) 



^(1) 



A7 



(0) 



Finally, the running coupling constant is given by 



47r 



1 



/5oln(/xVA' 



QCD. 



Aln[ln(/iVA^^^)] 
/32ln(/iVA|cD) 



+ ■ 



(44) 



(45) 



The above equation defines the continuum MS scale parameter Kqcd at the 
NLO. 

In the case of the A/ = 3/2 lattice operator, the coefficient and the 
anomalous dimension "-y^^^ become a vector due to the mixing of the bare 
operator with operators of different chirality, induced by the Wilson term 

. Operators of different chirality 



present in the lattice action [^, |26|, ^ 
correspond, in the language of refs. to the so-called "effervescent" 

operators. Denoting by T\ the zeroth-order Green function of the original 
operator and by Fq with i = 2, . . . 4 the zeroth-order Green functions of the 
"effervescent" operators [B^, , eq. (34) becomes 



F(H= E n + 



i=lA 



167r2 



5iA 



7 



(0) 



In 



7T 



C 



+ ( 



9lW 
167r2 



447(°H-2/?o + 7(°))ln^ {^y + lH% 



+ (-2/5o + 7(°))C^/^ + /3r (A 



3(0) 



dCt 
dX 



+ ■ ■ 



(46) 
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with \dCl/d\ ~ 

In the continuum, if we do not remove the W propagator and expand the 
four-fermion amphtude in powers of as{Mw), we do not need to introduce 
any regularization. In other words, the only divergences appearing in the 
four-fermion amplitude are those due to the renormalization of ag- Since the 
VT-propagator acts as an ultraviolet cut-off for the divergences of the Green 
function, we will use in the following the name Mv^-regularization, as much 
as we use the names dimensional or lattice regularization. We thus can write 



ro 1 



as{Mw) ( 7 



(0) 



47r 



In 



;7 



(0)/ 



Mw 



c 



(1) 



In 



w 



[1 
+ .. 



(47) 



where A = X{Mw). rdi/Mw) obeys the renormalization group equation 

(4 + Pio^s)^ + /?a(«.)A^ + 7(a.))r (t) = 0, 



with 



2 ' 



and 



das{t) 



1 dX{t) 



dt \{t) dt 

By solving the renormalization group equation for T{ij,/Mw), one finds 



(48) 

(49) 
(50) 

(51) 



+ 



asili) 



An 



Air 



■(- 



2Pi 



(52) 
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where = C""(^A(/i)^. A similar solution is found for r(/ia), where we 
encounter the complication due to the mixing with the effervescent operators 



r(H = E ro(- 



wa,(7r/a)N7(°V2/3or^ , a,(/i) 



j=l,4 



2^1 ^+ 4vr V 



(53) 



The comparison of eqs. ( |52D and (|53| ) allows us to write the effective Hamil- 
tonian in terms of the lattice bare operators 



Heff = E 0,{a)C,{Mwa), 

i=l,4 



(54) 



where 

C,{Mwa) 



c>i,i H 



+ 



Cis{T^/a) 



An 



2Pi 



An 



(55) 



Another convenient expression for Ci{Mwo,) is given by 

, /a,(Mty)N7(oV2/3or^ / , ^.(Mh/) - a,(7r/a) 7«/5o - 7^°^/5i 
Ci(Mp^a) = I — - , . ^i,i(H 



as(7r/a) 
47r ^ ' 



47r 



/ a,(MH/) x7'°V2/3o/ a,(MH.) - a,(7r/a) 7(^/^0 - 7WA 



An 



(56) 



Equation (56) has been obtained using the universality of the combinatior 



2/3o 



2A 







(57) 



^This relation is true only when we use the same coupling constant on the lattice and 
in the continuum. 
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Equation ( p5D is interpreted as follows: the continuum operator, defined in 
the Miy-regularization is matched into the lattice operator at the Mw scale 
through the factor (l + as{Mw)/ATc(^C'^ — C^)) and then evolved, according 
to eq. (|iT|), from the scale down to vr/a On the other hand, eq. 

(|56D is understood in the standard lattice language as follows: the lattice 
operator, is matched into the continuum operator at the scale vr/a through 
the factor — C^^ ^ and then evolved in the continuum, 

according to eq. (^81) , from the scale vr/a up to Myy p5|, Equations (^) 
and ( [56| ) apply with trivial changes to two-quark operators that renormalize 
multiplicatively on the lattice, or to operators that mix under renormalization 
also in the continuum. A few comments may be useful at this point: 

• The continuum, two-loop anomalous dimension in the Mvi/-regularization 
7^^) can be obtained from the calculation done in the NDR, HV or 
DRED regularizations which we denote as 

(^Ms)' using 

the relation 

• In eq. (|56[), the coefficients C/" can be eliminated by defining a new 
lattice operator 

0\a) = Oi(a) - E CfO.(a). (59) 

i=l,4 

This is equivalent to the regularization-independent definition of the 
renormalized operators discussed in refs. and in sec. ^ This 

definition is such that r(/i = vr/a) = Fq []. We stress again that this 
procedure requires the knowledge of the external states for which we 
have computed and it is in general gauge- dependent. 

• We have decided to expand the lattice Green function in ^^(vr/a). How- 
ever, we have the freedom to expand in ^^(l/a) or any other scale we 



^ This is not in contrast with the condition <C [tt / of' that we have to impose in order 
to avoid discretization errors. The renormahzation condition is simply a consequence of eq. 
(34), which is derived by expanding the Green function in Q;s(7r/a), instead of g\{a)/AT:, 
with < (tt/o)^. 
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like, since the change will be completely compensated at the NLO For 
the same reason we can expand the Green functions in a different cou- 
phng constant, for example the "boosted" coupling defined in refs. 
]55| . There it was argued that the series in may minimize 0{al) 



NNLO corrections. If we expand in , we have accordingly to reorga- 
nize the expression used for the coefficient function. It is wise however 
to check the stability of the physical amplitudes under a change in the 
scale used for as and assume this as a theoretical uncertainty. 



From eq. (|5q) it is possible to give the expression of the operator from 
which we can immediately derive the renormalization group-invariant 
Bk parameter. This is obtained by using the operator 

Orgi = C{a)0{a), (60) 

with 

0(a) = 0,(.)(l + ^CO) - -J^Y^C-0,(a) (61) 

or, in the MS scheme. 



0(a)^ = 0,(a)(l + '^M^C^^) - V C^OM) (62) 

V 471- ^ A-K 



1=1,4 



and 

r( \ ( ( I ^r^*°'/'*h «s(7r/a) ^ (7^^)/9o-7^°^A M .... 
C[a) = [as{7r/a)) [1 jj ; (63) 



C^^ can be found in refs. g, |§ . 

Contrary to what is stated in ref. |Q, it is not necessary to first 
match the lattice operators to the corresponding operators in some con- 
tinuum renormalization scheme at a low scale /xq, with the W integrated 
out, and then evolve to Mw- The reason is that, as demonstrated 
above, the theory with the W does not require any renormalization 



® This is true provided that the scale that we choose is not so different from 1/a as to 
generate new large logarithms. 
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when expanded in as{Mw)- Thus we can avoid the introduction of /io 
and relate directly the bare lattice operators to the full theory. The 
standard approach to evaluate the theoretical uncertainty, coming from 
the matching of the lattice to the continuum, is obtained by varying 
fiQ and the effective ("boosted") lattice coupling constant, used in the 
perturbative expansion. We have shown that, indeed, the uncertainty 
depends only on the effective coupling, since in our formulae /xq never 
appears. One can use the operator 0'(a) of eq. (|59D, computed using 
the non-perturbative procedure of ref. [^. In this case, the matching 
can be entirely done in the continuum, so that the large corrections, 
present in higher orders in the lattice coupling, are eliminated. 



6 The RI renormalization scheme 

The details of the renormalization of the operators and of the calculation of 
the Wilson coefficients in HV and NDR have been extensively discussed in 
the literature 0-0. Here, we give the main formulae, necessary to renormal- 
ize the operators in the RI renormalization scheme. Using these formulae, 
the reader can obtain the renormalized operators, starting from any regular- 
ized version of QCD. The matrix elements of the RI renormalized operators 
can be combined directly with the coefficients given in tables 4 and 5. 

At the end of this section, we will give the practical rules to be followed, 
in order to obtain the physical amplitudes. 



The construction of the renormalized operators appearing in (0) is com- 
plicated by their mixing and by the presence of the so-called "effervescent" 
operators. To simplify the discussion, we start by considering the A J = 3/2 
(AS* = 2) operator, which, in the absence of "effervescent" operators, renor- 
malizes multiplicatively. The four-fermion A/ = 3/2 operator mixes at one 
loop with other ("effervescent") operators, due to the artefacts of the reg- 
ularization. This problem is common to the continuum and to the lattice 
cases^ The first step is to introduce counter-terms to eliminate this artificial 
mixing. This can be achieved in many different ways. On the lattice, for 

^"Effervescent" operators appear in all regularizations in the presence of 75. They are 
also present using the procedure of ref. ||57| , which avoids to define 75 . 
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Feynman RI 






Landau RI 












= 1.5 GcV 








(-2.72 ± 0.29 ± 0.00 ± 0.06) 


• 10 


-i 


(-3.33 ± 0.44 ± 0.00 ± 0.04) 


■ 10 


-1 


C2 


(108.70 ± 0.73 ± 0.00 ± 0.25) 


• 10 


-2 


(11.59 ±0.24 ±0.00 ±0.06) 


10" 


-1 


Ca 


(2.05 ±0.27 ±0.00 ±0.11) ■ 


10- 


2 


(2.35 ±0.38 ±0.00 ±0.07) • 


10" 


2 


Ci 


(-5.60 ± 0.68 ± 0.01 ± 0.10) 


• 10 


-2 


(-5.93 ±0.81 ±0.01 ±0.05) 


■ 10 




Cs 


(10.00 ± 0.91 ± 0.01 ± 0.50) 


10" 


-3 


(11.54 ±0.40 ±0.01 ±0.42) 


10" 




Ce 


(-1.17 ± 0.30 ± 0.01 ± 0.07) 


• 10 


-1 


(-1.12 ± 0.28 ± 0.01 ± 0.06) 


• 10 




C7 


(0.02 ± 0.07 ± 0.18 ± 0.02) ■ 


10" 


3 


(0.01 ±0.07 ±0.19 ±0.02) • 


10" 


3 


Cs 


(1.68 ± 0.48 ± 0.30 ± 0.08) • 


10" 


3 


(1.56 ± 0.44 ± 0.28 ± 0.09) • 


10" 


3 


C9 


(-7.43 ± 0.21 ± 0.75 ± 0.14) 


■ 10 


-3 


(-7.54 ± 0.19 ± 0.76 ± 0.13) 


• 10 


—3 






M 


= 2. 


GeV 






Ci 


(-2.27 ± 0.22 ± 0.00 ± 0.04) 


■ 10 


-i 


(-2.75 ± 0.32 ± 0.00 ± 0.03) 


• 10 


— 1 


C2 


(106.88 ± 0.60 ± 0.00 ± 0.21) 


■ 10 


-2 


(11.27 ± 0.17 ± 0.00 ± 0.04) 


10" 


- 1 


C3 


(1.73 ± 0.20 ± 0.01 ± 0.07) • 


10~ 


2 


(1.93 ±0.26 ±0.01 ±0.04) • 


10- 


2 


04 


( 4 -t- n ^^9 -k n nn -k n nf^^ 

1^ 'i.o'i HI u.oz zn. u.uu zn. u.uuj 


. 1 n 

-LU 


-2 


(-5.07 ±0.59 ±0.01 ±0.03) 


-LU 


-2 


C5 


(^lU.oi ± U.lo it U.Uz it U.oiJ 


iU 


-3 


(11.43 ±0.26 ±0.02 ±0.25) 


1 

lU 


-3 


^6 


— u.yu ^ u.xo ^ u.uu ^ u.ui^ 


■ 10 


-1 


(-0.86 ±0.16 ±0.00 ±0.04) 


■ 10 


-1 




/■fi nn -U n n/i -U n i o -U n no^ 
I^U.Uu it U.u4 it U.iy it \j.\)jL) ■ 


1 n — 
iU 


3 


(-0.01 ±0.04 ±0.20 ±0.02) 


1 n 
- iU 


-3 




( ^ ofi_l_no'7_l_^0'3J_nnc^ 
[i.Zd it U.Zi it U.Zo it U.UoJ ■ 


in — 
iU 


3 


(1.18 ±0.25 ±0.21 ± 0.09) • 


1 n — 
iU 


3 


( ■!) 


/ 7 vi _i_ n "1 -U n TQ -U n 1 fi^i 

— (./ i zt U.IO zt U./o zt U.iO) 


■ lu 


-3 


(-7.81 ±0.14 ±0.79 ±0.16) 


1 n 


-3 






/' 


= :-i.O Gq\' 






Ci 


(-1.74 ± 0.15 ± 0.00 ± 0.03) 


• 10 


-i 


(-2.11 ± 0.21 ± 0.00 ± 0.02) 


• 10 


— i 


C2 


(104.70 ± 0.38 ± 0.00 ± 0.18) 


■ 10 


-2 


(10.92 ± 0.10 ± 0.00 ± 0.03) 


10" 


- 1 


C3 


(1.37 ± 0.13 ± 0.01 ± 0.05) • 


10~ 


2 


(1.50 ±0.16 ±0.01 ±0.03) • 


10- 


2 


C4 


(-3.97 ± 0.37 ± 0.00 ± 0.04) 


• 10 


-2 


(-4.12 ± 0.41 ± 0.00 ± 0.02) 


■ 10 


-2 


C5 


(9.78 ±0.35 ±0.03 ±0.18) ■ 


10- 


3 


(10.53 ± 0.52 ± 0.03 ± 0.13) 


10" 


-3 


Ce 


(-6.49 ±0.96 ±0.03 ±0.25) 


■ 10 


-2 


(-6.29 ±0.91 ±0.03 ±0.21) 


■ 10 


-2 


C7 


(0.00 ±0.02 ±0.19 ±0.02) ■ 


10- 


3 


(0.01 ±0.02 ±0.20 ± 0.02) ■ 


10- 


3 


Cs 


(0.90 ±0.15 ±0.16 ±0.08) ■ 


10" 


3 


(0.84 ±0.14 ±0.15 ±0.08) • 


10" 


3 


C9 


(-7.93 ± 0.11 ± 0.81 ± 0.18) 


■ 10 


-3 


(-8.01 ± 0.10 ± 0.82 ± 0.18) 


■ 10 


-3 








= 4.0 GeV 






Ci 


(-1.41 ± 0.12 ± 0.00 ± 0.02) 


■ 10 


-i 


(-1.73 ± 0.16 ± 0.00 ± 0.01) 


• 10 


-i 


C2 


(103.40 ± 0.26 ± 0.00 ± 0.17) 


• 10 


-2 


(107.32 ± 0.74 ± 0.00 ± 0.29) 


• 10 


-2 


C3 


(11.80 ± 0.98 ± 0.14 ± 0.35) 


10" 


-3 


(1.27 ±0.12 ±0.01 ±0.02) • 


10" 


2 


C4 


(-3.46 ± 0.29 ± 0.00 ± 0.03) 


• 10 


-2 


(-3.57 ± 0.32 ± 0.00 ± 0.01) 


• 10 


-2 


Cs 


(9.12 ±0.43 ±0.03 ±0.13) • 


10" 


3 


(9.70 ± 0.54 ± 0.03 ± 0.09) • 


10" 


3 


Cs 


(-5.28 ± 0.67 ± 0.03 ± 0.18) 


■ 10 


-2 


(-5.14 ± 0.64 ± 0.03 ± 0.15) 


■ 10 


-2 


C7 


(0.03 ±0.01 ±0.19 ±0.02) • 


10- 


3 


(0.04 ±0.01 ±0.20 ±0.02) ■ 


10- 


3 


Cs 


(0.73 ±0.10 ±0.14 ±0.07) ■ 


10- 


3 


(0.69 ±0.10 ±0.12 ±0.07) ■ 


10- 


3 


Cg 


(-8.11 ±0.08 ±0.83 ±0.19) 


■ 10 


-3 


(-8.18 ± 0.08 ± 0.84 ± 0.19) 


■ 10 


-3 



Table 4: Wilson coefficients of the effective AS = 1 Hamiltonian in RI for 
A* = and 1. The first error comes from the uncertainty on Aqcd, the 
second one from that on the top mass. The third error is an estimate of 
the uncertainty coming from the use of different formulae for the coefficient 
functions, which are equivalent at NLO order. 
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Feynman RI 




Landau RI 








^l = 4.5 GeV 


Ci 


(-1.29 ±0.10 ±0.00 ±0.02) 


■10-^ 


(-1.58 ±0.14 ±0.00 ±0.01) 


■ 10 


-i 




(102.92 ± 0.21 ± 0.00 ± 0.17) 


■ 10-2 


(106.63 ± 0.64 ± 0.00 ± 0.27) 


■ 10 


-2 


Ca 


(12.38 ±0.97 ±0.32 ±0.29) 


10-3 


(1.35 ±0.12 ±0.03 ±0.02) 


10- 


2 


C4 


(-3.40 ± 0.28 ± 0.02 ± 0.03) 


■ 10-2 


(-3.53 ±0.30 ±0.02 ±0.01) 


■ 10 


-2 


Cs 


(8.78 ±0.42 ±0.03 ±0.11) ■ 


10-:^ 


(9.31 ±0.52 ±0.03 ±0.06) 


10- 


3 


Ce 


(-4.88 ±0.59 ±0.03 ±0.16) 


■ 10-2 


(-4.75 ±0.56 ±0.03 ±0.13) 


■ 10 


-2 


C7 


(0.05 ±0.01 ±0.19 ±0.02) • 


10-3 


(0.05 ± 0.01 ± 0.20 ± 0.02) 


10- 


3 


Cs 


(0.68 ± 0.09 ± 0.12 ± 0.07) • 


10-3 


(0.63 ±0.08 ±0.12 ±0.07) 


10- 


3 


C9 


(-0.95 ± 0.00 ± 0.10 ± 0.02) 


• 10-2 


(-0.98 ± 0.01 ± 0.10 ± 0.02) 


• 10 


-2 


Cio 


(1.34 ±0.10 ±0.17 ±0.03) • 


10-3 


(1.60 ±0.13 ±0.20 ±0.03) 


10" 


3 




M = 5.0 GeV 


Ci 


(-11.85 ± 0.94 ± 0.00 ± 0.17) 


■ 10-2 


(-1.46 ±0.13 ±0.00 ±0.01) 


■ 10 


-1 


C2 


(102.50 ± 0.17 ± 0.00 ± 0.16) 


■ 10-2 


(106.05 ±0.57 ±0.00 ±0.26) 


■ 10 


-2 


C3 


(11.75 ±0.89 ±0.32 ±0.27) 


10-3 


(1.27 ±0.11 ±0.03 ±0.02) 


10- 


2 


C4 


(-3.25 ± 0.26 ± 0.02 ± 0.02) 


■ 10-2 


(-3.36 ±0.28 ±0.02 ±0.01) 


■ 10 


-2 


Cs 


(8.53 ±0.43 ±0.03 ±0.10) ■ 


10-3 


(9.03 ±0.52 ±0.03 ±0.05) 


10- 


3 


Ce 


(-4.57 ±0.53 ±0.03 ±0.14) 


■ 10-2 


(-4.45 ±0.50 ±0.03 ±0.12) 


■ 10 


-2 


C7 


(0.06 ± 0.01 ± 0.20 ±0.02) • 


10-3 


(0.07 ±0.01 ±0.20 ±0.02) 


10- 


3 


Cg 


(0.63 ±0.08 ±0.12 ±0.07) • 


10-3 


(0.59 ±0.07 ±0.11 ±0.07) 


10- 


3 


O) 


(-0.91 ±0.00 ±0.10 ±0.02) 


• 10-2 


(-0.97 ± 0.01 ± 0.10 ± 0.02) 


• 10 


-2 


Ciu 


(1.27 ± D.OI) ± O.K) ± 0.0:-!) • 


10 


(1.52 ± 0.12 ± 0.19 = 0.0:-!) 


10 






fi = 10.0 GcV 


Ci 


(-5.68 ± 0.41 ± 0.00 ± 0.10) 


■ 10-2 


(-7.78 ± 0.61 ± 0.00 ± 0.07) 


■ 10 


-2 


C2 


(100.19 ± 0.02 ± 0.00 ± 0.13) 


■ 10-2 


(102.95 ± 0.23 ± 0.00 ± 0.19) 


■ 10 


-2 


C3 


(8.28 ± 0.51 ± 0.30 ±0.16) • 


10-3 


(8.84 ±0.60 ±0.32 ±0.09) 


10- 


3 


C4 


(-2.35 ± 0.17 ± 0.02 ± 0.01) 


■ 10-2 


(-2.41 ±0.18 ±0.02 ±0.00) 


■ 10 


-2 


C5 


(6.77 ±0.39 ±0.03 ±0.05) ■ 


10-3 


(7.05 ±0.43 ±0.04 ±0.02) 


10- 


3 


Ce 


(-2.97 ± 0.27 ± 0.02 ± 0.08) 


■ 10-2 


(-2.91 ±0.26 ±0.02 ±0.07) 


■ 10 


-2 


C7 


(0.20 ±0.00 ±0.20 ±0.02) ■ 


10-3 


(0.21 ±0.00 ±0.21 ±0.02) 


10- 


3 


C8 


(4.17 ±0.41 ±0.77 ±0.55) ■ 


10-4 


(3.87 ±0.38 ±0.70 ±0.55) 


10- 


4 


C9 


(-9.11 ±0.01 ±0.98 ±0.17) 


• 10-3 


(-0.94 ± 0.00 ± 0.10 ± 0.02) 


• 10 


-2 


Cio 


(7.67 ± 0.47 ± 0.99 ± 0.45) • 


lQ-4 


(0.95 ±0.07 ±0.12 ±0.05) 


10" 


3 



Table 5: Same as in table 4 o,t different values of the scale jjL. 
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example, it is obtained by using eq. (^9]). In dimensional regularizations, the 
elimination of the unwanted operators is ensured by minimal subtraction of 
the pole term: the string of gamma matrices, proportional to 1/e at one loop, 
cancels all "effervescent" contributions 0-0. There is some ambiguity in 
the subtraction. In eq. (|59[) , one can, for example, subtract all the terms 
but CfOi(a) or, in dimensional regularization, one can subtract any further 
piece of 0(e°) proportional to the original four-fermion operator. The impor- 
tant point is that, after this first step, the subtracted operator Osub is only 
proportional to the original one. Notice that, at one loop, the coefficients 
of the "effervescent" operators are independent of the external states and 
consequently they are gauge- independent, as explicit one- loop calculations 
show. The independence from the external states ensures the possibility of 
defining renormalized operators with definite chiral properties, by combining 
bare and "effervescent" operators in a gauge-invariant way. At this point, 
we impose the overall renormalization condition to the subtracted four-quark 
Green function To^^^ip), by taking = fi'^ for all the external legs[^ Let us 
define 



= 1 + ^ (-7^°) HpVf^') + r""'^'''''' + (1 - A)rA) , (64) 

where P is a suitable projector operator and A denotes the gauge 

parameter defined in the previous section; Aq ^(p) includes the renormal- 
ization of the external quark legs. In eq. (|64[), r^^'^^^ depends on the 
regularization in which A^ (p) is calculated, whereas r\ is regularization- 
and gauge- independent. On the lattice, eq. (0) has a similar expression with 
ln(p^//i^) — > ln(pa/7r)^ and r^^'^^^ —>■ r^. The renormalization condition 

{Zg\f^, Xny' X A^^^^(p)|,.=^.,,=,. = 1 (65) 
defines the renormalized operator 

0^^(/i, A*) = {Zg'ifi, X*)y' Osub. (66) 

^''To avoid unnecessary complications, we take all the external momenta equal to fP, 
i.e. the scale that multiplies the coupling constant in dimensional schemes. Of course, one 
can choose arbitrary external momenta, provided they make the Green functions infrared- 
finite. In the present discussion, we assume that we work in the Euclidean space, in order 
to make contact with the previous section. 
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Then, one has 



Ao-(P) = 1 + (-7^°^ ln(pV/i') + (A* - X)rx) , (67) 

where A* is the gauge introduced in eq. (^); A is instead the gauge in which 
the Green function, containing the insertion of 0^\fi, A*), is computed. 

If we call r^^ the finite part of the matrix element of O^^ , from eq. (|571) 



we find 



Equation (|68|) has a very simple interpretation: the first term, —fHV,NDR^ 
removes the regularization dependence; the second one, ta, introduces a de- 
pendence on A*; rx depends on the external states but not on the regular- 
ization. The ambiguity of the first step, i.e. in the subtraction of the pole 
term, discussed above, is removed by the term —j-^^^^^^^ which depends on 
the subtraction procedure. 

Almost all the perturbative lattice calculations have been done in the 
Feynman gauge, while the Landau gauge is the most convenient one to im- 
plement the non-perturbative method of ref. For this reason, we give 
the Wilson coefficients for A* = 0, 1. 

In tables 4-7, the coefficients are given at several values of /x, in the 
RI (A* = 0, 1), NDR and HV renormalization schemes. We have extended 
the range of /x, at which the coefficients are computed, up to 10 GeV, so 
that they can also be used for 5-decays. We present separately the error 
coming from the uncertainty on Aqcd, on the value of the top quark mass. 
In computing the coefficients, several expressions, which are equivalent at 
the order at which we are working, can be used. For this reason, one may 
find smallish differences (of 0(a^)) between our results and those presented 
in ref. 0. The third error in tables 4 and 5 is an estimate of these O(a^) 
effects, that have been evaluated by using for the Wilson coefficients different 
expressions which are equivalent at the NLO. 

We summarize the RI renormalization prescription 

1. Define a subtracted operator by removing the pole terms and the mix- 
ing with "effervescent" operators. Any further finite subtraction pro- 
portional to the original operator is immaterial. As explained in the 
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LO 


NLO HV 


NLO NDR 




/^ = 1.5 GeV 


C2 
C3 

Ca 
Cs 
Ce 
Ct 
Cs 
Co 


(-4.22 ± 0.65 ± 0.00) x 10"^ 
(11.62 ± 0.38 ± 0.00) X 10-1 

(1.99 ±0.35 ±0.00) X 10-2 
(-4.16 ± 0.56 ± 0.02) X 10-^ 

(1.19 ±0.12 ±0.00) X 10-2 
(-0.66 ± 0.13 ± 0.00) X 10-1 

(0.16 ±0.04 ±0.19) X IQ-s 

(0.63 ±0.14 ±0.16) X 10-^ 
(-6.77 ±0.27 ±0.71) x 10-^ 


(-3.91 ± 0.51 ± 0.00) X 10-1 
(106.13 ± 0.82 ± 0.00) X 10-2 

(2.17 ±0.41 ±0.00) X 10-2 
(-4.51 ± 0.60 ± 0.01) X 10-2 

(1.37 ±0.15 ±0.00) X 10-2 
(-0.63 ± 0.11 ± 0.00) X 10-1 
(-0.04 ± 0.00 ± 0.17) X 10-3 

(0.97 ±0.19 ±0.15) X 10-3 
(-6.32 ± 0.37 ± 0.64) x 10-3 


(-3.80 ± 0.55 ± 0.00) X lO-i 
(11.95 ± 0.35 ± 0.00) X 10-1 
(2.60 ±0.52 ±0.00) X 10-2 
(-0.63 ± 0.11 ± 0.00) X 10-1 
(10.52 ± 0.61 ± 0.01) X 10-3 
(-0.93 ± 0.21 ± 0.00) X 10-1 
(0.02 ±0.06 ±0.20) X 10-3 
(1.06 ±0.26 ±0.19) X 10-3 
(-7.24 ± 0.19 ± 0.73) x 10-3 




M = 2 GeV 


Ci 

C2 

C3 
C4 
Cs 
Cs 
C7 
Cs 
C9 


(-3.47 ± 0.44 ± 0.00) x 10"^ 
(11.16 ±0.23 ±0.00) X 10-1 
(1.59 ±0.23 ±0.00) X 10-2 
(-3.50 ± 0.40 ± 0.01) X 10-2 
(10.40 ± 0.94 ± 0.04) X 10-^ 
(-5.23 ± 0.80 ± 0.03) x IO-2 
(0.18 ±0.02 ±0.19) X 10-3 
(0.50 ±0.08 ±0.12) X 10-=* 
(-7.03 ± 0.20 ± 0.74) x 10"'^ 


(-3.29 ± 0.37 ± 0.00) x IQ-i 
(104.13 ± 0.54 ± 0.00) X 10-2 

(1.73 ±0.26 ±0.00) X 10-2 
(-3.82 ± 0.44 ± 0.01) X 10-2 

(1.20 ±0.11 ±0.00) X 10-2 
(-5.08 ±0.72 ±0.03) x IO-2 

(0.01 ±0.00 ±0.18) X 10-3 

(0.77 ±0.12 ±0.12) X 10-3 
(-6.71 ±0.27 ±0.68) x 10-3 


(-3.13 ± 0.39 ± 0.00) X 10-1 
(11.54 ±0.23 ±0.00) X 10-1 
(2.07 ±0.33 ±0.00) X IO-2 
(-5.19 ± 0.71 ± 0.01) X 10-2 
(10.54 ± 0.16 ± 0.02) X 10-3 
(-0.72 ± 0.13 ± 0.00) X 10-1 
(0.01 ±0.04 ±0.20) X 10-3 
(0.81 ±0.16 ±0.14) X 10-3 
(-7.49 ± 0.15 ± 0.75) x 10-3 




^ = 3 GcV 


Ci 
C2 
C3 
C4 
Cs 
Cs 
C7 
Cs 
C9 


(-2.68 ±0.28 ±0.00) x 10"^ 
(10.71 ±0.12 ±0.00) X 10-1 

(1.20 ±0.14 ±0.01) X 10^2 
(-2.78 ± 0.26 ± 0.01) X 10^2 

(8.60 ± 0.67 ± 0.04) x 10-^ 
(-3.89 ± 0.47 ± 0.03) x IO-2 

(0.22 ±0.01 ±0.20) X 10-3 

(3.73 ± 0.50 ± 0.92) x lO-'' 
(-7.29 ± 0.14 ± 0.77) x 10-^ 


(-2.59 ± 0.25 ± 0.00) x IQ-i 
(101.88 ± 0.24 ± 0.00) X IO-2 

(1.29 ±0.16 ±0.00) X 10-2 
(-3.06 ±0.29 ±0.01) X 10-2 

(9.99 ± 0.84 ± 0.04) x 10-3 
(-3.84 ± 0.44 ± 0.03) x IO-2 

(0.08 ± 0.00 ± 0.19) X 10-3 

(5.84 ±0.74 ±0.91) x 10-"* 
(-7.04 ± 0.18 ± 0.73) X 10-3 


(-2.41 ±0.25 ±0.00) X 10-1 
(11.12 ±0.14 ±0.00) X 10-1 

(1.56 ±0.19 ±0.01) X 10-2 
(-4.10 ± 0.46 ± 0.00) X 10-2 

(9.75 ± 0.46 ± 0.03) x 10-3 
(-5.33 ± 0.73 ± 0.03) x IO-2 

(0.04 ±0.02 ±0.20) X 10-3 

(0.59 ±0.09 ±0.10) X 10-3 
(-7.66 ± 0.11 ± 0.78) X 10-3 




ft = 4.0 GcV 


Ci 
C2 
C3 
C4 
Cs 
Cs 
C7 
Cs 
C9 


(-2.23 ±0.21 ±0.00) x 10-1 
(104.69 ± 0.80 ± 0.00) x 10^2 

(9.91 ±0.98 ±0.09) x lO^^ 
(-2.36 ± 0.20 ± 0.00) x 10^2 

(7.48 ± 0.54 ± 0.04) x 10-3 
(-3.18 ± 0.33 ± 0.02) X 10-2 

(0.27 ±0.01 ±0.20) X 10-3 

(3.10 ± 0.36 ± 0.75) X 10-* 
(-7.44 ± 0.11 ± 0.79) X 10-3 


(-2.18 ±0.19 ±0.00) X 10-1 
(1005.88 ± 0.92 ± 0.00) x 10-3 

(1.06 ±0.11 ±0.01) X 10-2 
(-2.61 ±0.23 ±0.00) X 10-2 

(8.75 ± 0.69 ± 0.04) x 10-3 
(-3.18 ± 0.32 ± 0.02) X 10-2 

(0.15 ±0.00 ±0.19) X 10-3 

(4.87 ± 0.55 ± 0.75) x 10-* 
(-7.27 ± 0.15 ± 0.76) x 10-3 


(-1.99 ±0.19 ±0.00) X 10-1 
(109.00 ± 0.99 ± 0.00) X IO-2 

(1.29 ±0.14 ±0.01) X 10-2 
(-3.50 ± 0.35 ± 0.00) X IO-2 

(8.97 ± 0.50 ± 0.03) x 10-3 
(-4.38 ± 0.51 ± 0.02) X 10-2 

(0.08 ±0.01 ±0.20) X 10-3 

(4.76 ± 0.60 ± 0.84) x 10-* 
(-7.83 ± 0.09 ± 0.80) x 10-3 



Table 6: Wilson coefficients for the HV and NDR MS schemes. The first 
error comes from the uncertainty on Aqcd, the second one from that on the 
top mass. 
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LO 


NLO HV 


NLO NDR 




fj, = 4.5 GcV 


Ca 
Ca 
Ci 
Cs 
Ce 
Ct 
Cs 
C9 
Cio 


(-2.07 ±0.19 ±0.00) X 10-^ 
(103.85 ±0.67 ±0.00) x 10"^ 

(1.13 ± 0.11 ± 0.03) X 10-2 
(-2.42 ± 0.20 ± 0.02) X 10"^ 

(7.06 ± 0.49 ± 0.04) x lO^^ 
(-2.94 ± 0.29 ± 0.02) X lO'^ 

(0.29 ± 0.01 ± 0.20) X 10-3 

(2.88 ± 0.32 ± 0.69) x 10-'' 
(-0.96 ± 0.01 ± 0.10) X 10-2 

(2.12 ± 0.19 ± 0.24) X 10-3 


(-2.03 ±0.18 ±0.00) X lO--' 
(1001.15 ± 0.43 ± 0.00) X 10-3 

(1.19 ±0.12 ±0.03) X 10-2 
(-2.67 ± 0.23 ± 0.02) x IO-2 

(8.36 ±0.66 ±0.04) x 10-3 
(-2.97 ± 0.29 ± 0.02) x IO-2 

(0.17 ±0.00 ±0.19) X 10-3 

(4.53 ± 0.49 ± 0.70) x 10-" 
(-0.93 ± 0.00 ± 0.10) x 10-2 

(1.95 ±0.17 ±0.24) X 10-3 


(-1.84 ±0.17 ±0.00) X 10-^ 
(108.20 ± 0.87 ± 0.00) x IO-2 

(1.40 ±0.14 ±0.03) X 10-2 
(-3.51 ± 0.34 ± 0.02) X IO-2 

(8.75 ± 0.52 ± 0.04) x 10-3 
(-4.09 ± 0.46 ± 0.02) x IO-2 

(0.09 ±0.01 ±0.20) X 10-3 

(4.39 ± 0.52 ± 0.77) x 10-* 
(-0.97 ± 0.01 ± 0.10) X 10-2 

(1.86 ±0.16 ±0.23) X 10-3 




^ = 5 GeV 


Ci 

Ca 
C3 
C4 
Cs 
Ce 
C7 
Cs 
C, 
Cio 


(-1.93 ± 0.17 ± 0.00) X 10"^ 
(103.17 ± 0..56 ± 0.00) X 10"^ 
(10.60 ± 0.95 ± 0.34) x 10^3 
(-2.29 ± 0.18 ± 0.02) X 10"^ 

(6.72 ± 0.46 ± 0.04) x 10^3 
(-2.74 ± 0.26 ± 0.02) X 10"^ 

(0.31 ± 0.01 ± 0.20) X 10-3 

(2.70 ±0.29 ±0.65) x 10-* 
f-0.95 ± 0.01 ± 0.10) X 10"'-' 

(2.01) ± 0.18 ± 0.2:-!) X 10 


(-1.90 ±0.16 ±0.00) X 10-^ 
(9970.37 ± 0.31 ± 0.00) x 10-* 

(1.11 ±0.10 ±0.03) X 10-2 
(-2.53 ±0.21 ±0.02) X IO-2 

(7.97 ±0.62 ±0.04) x 10-3 
(-2.78 ±0.26 ±0.02) x IO-2 

(0.20 ±0.00 ±0.19) X 10-3 

(4.25 ±0.44 ±0.65) x 10-* 
(-9.23 ± 0.02 ± 1.00) x 10-3 

(1.87 = 0.15 = 0.22) X 10 


(-1.71 ±0.15 ±0.00) X 10-^ 
(107.54 ± 0.77 ± 0.00) X IO-2 

(1.32 ±0.13 ±0.03) X 10-2 
(-3.32 ± 0.31 ± 0.02) X IO-2 

(8.45 ±0.52 ±0.04) x 10-3 
(-3.83 ± 0.41 ± 0.02) X IO-2 

(0.11 ±0.01 ±0.20) X 10-3 

(4.09 ± 0.46 ± 0.71) X 10-* 
(-0.96 ±0.01 ±0.10) X 10-2 

(1.77 ±0.11 ±0.22) X 10 3 




^ = 10 GcV 


Ci 
Ca 
C3 
C4 
Cs 
Ce 
C7 
Cs 
C9 
Cio 


(-11.76 ± 0.90 ± 0.00) X 10-^ 
(99.56 ±0.11 ±0.00) X 10-2 

(6.86 ± 0.49 ± 0.32) x 10^3 
(-1.52 ± 0.11 ± 0.02) X 10^2 

(4.67 ±0.30 ±0.04) x 10^3 
(-1.71 ± 0.14 ± 0.02) X 10^2 

(0.46 ± 0.00 ± 0.20) X 10^3 

(1.76 ± 0.15 ± 0.40) X 10-"' 
(-9.06 ± 0.04 ± 1.00) X 10-3 

(1.33 ±0.10 ±0.15) X 10-3 


(-11.74 ± 0.88 ± 0.00) X 10-2 
(97.47 ±0.18 ±0.00) x IO-2 
(7.12 ±0.53 ±0.31) X 10-3 

(-1.69 ±0.13 ±0.02) X 10-2 
(5.56 ± 0.40 ± 0.04) x 10-3 

(-1.75 ± 0.14 ± 0.02) X 10-2 
(0.39 ±0.01 ±0.20) X 10-3 
(2.74 ±0.23 ±0.41) X 10-* 

(-8.87 ± 0.02 ± 0.98) x 10-3 
(1.24 ±0.09 ±0.15) X 10-3 


(-9.73 ± 0.77 ± 0.00) x 10-^ 
(104.01 ± 0.34 ± 0.00) X 10-2 

(8.74 ±0.67 ±0.32) x 10-3 
(-2.29 ± 0.18 ± 0.02) X IO-2 

(6.42 ±0.41 ±0.04) X 10-3 
(-2.48 ± 0.22 ± 0.02) X 10-2 

(0.27 ±0.00 ±0.20) X 10-3 

(2.57 ±0.23 ±0.44) x 10-* 
(-9.19 ± 0.04 ± 1.00) X 10-3 

(1.17 ±0.08 ±0.15) X 10-3 



Table 7: Same as in table 6 hut at different scales ijl. 
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p 



Figure 5: Penguin diagram. We show explicitly the external momenta used 
in the RI renormalization scheme. 



previous section, by working with bare lattice operators, one only needs 
to eliminate the "effervescent" terms. 



2. Impose the renormalization condition ([65| ) to the subtracted operator 
in a fixed gauge A*. 

3. Take the coefficients from table 4 or 5 (when A* = 0, 1). Notice that 
the non-perturbative method cannot be applied in the Coulomb 
gauge. For other gauges some more analytic work is needed. 



In the more complicated case of operator mixing, i.e. mixing between 
operators belonging to the physical basis, e.g. (|15]), one may proceed in 
close analogy, by imposing that each renormalized operator of the physical 
basis has non-zero projection only onto itself. For penguin diagrams, we 
cannot take all external legs with the same momentum. We then choose the 
convenient configuration of momenta given in fig. ^. For completeness the 
matrices Ar^^'^^^, for A* = 0, 1 are given in table 8. 
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Table 8: 



The matrices Ar^X=o 1 ^'^^ ^'"^■=^1 necessary to relate the Wilson 



coefficients in the HV and NDR MS schemes to those computed in the RI 
scheme, in the Landau (X* = 0) and Feynman (X* = 1) gauges, are given. 
{i,j) refers to the matrix element. The matrix elements which are zero are 
not reported. ( Continue ) 
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^' A*=0 


^' A*=l 


^' A'=0 


^' A*=l 


(10,3) 


bu 1 bd 


bu 1 5d 
27 54 


5it 1 5d 


5u 1 5d 


27 ^ 54 


27 54 


27 ~^ 54 


(10,4) 


bu 5a 


bu bd 


5u 5a 


5ii 5a 


9 18 


9 18 


9 18 


9 18 


(10,5) 


bu I bd 


bu 1 5a 
27 54 


bu I 5a 


bu I 5a 


27 54 


27 54 


27 54 


(10,6) 


bu bd 


bu bd 


bu bd 


bu bd 


9 18 


,,9 18 


9 18 


,-9 18 


(10,9) 


5 - 121n2 


^ - 161n 2 

_iz 4- 16 In 2 


7- 121n2 


^ - 161n2 


(10,10) 


-7 + 41n2 


+ 41n2 


-| + fin 2 



Table 8: (Continued) 



7 5-parameters of the relevant operators 



Since our last study of e'/e [0, no significant progress, in lattice calculations 
of hadronic matrix elements, has been made. We want, however, to discuss 
a few points not clarified in our previous papers or raised by other authors. 



1. In ref. |^, for the kaon i?-parameter in NDR, they found BxifJ' = 2 
GeV) = 0.616 ± 0.020 ± 0.017. Using this value and the LO con- 
version factor, they derive the renormalization group-invariant Bk = 
0.825 ± 0.027 ± 0.023. For consistency with the present application, 
since corrections of 0{as) were included in the lattice calculation, we 
have to use the NLO expression (^) instead. In this case, one obtains 
Bk = 0.75 as central value. We have thus used B^ = 0.75, assuming an 
error of about 0.15, because of uncertainties coming from the quenched 
approximation, the scale in the coupling constant, Aqcd, etc. 

2. In the absence of electromagnetic corrections, it is possible to demon- 
strate (non-perturbatively) that 

B',^\f^) = BM (69) 

in any regularization scheme and at any scale fi The proof con- 
sists in showing that, by assuming SU{3) flavour symmetry, the ma- 
trix elements (7r+7r°|Q9|i^+)A7=3/2 and {n+7T^\Q+\K+)Ai=3/2 (Q+ = 
l/2{Qi+Q2)) correspond to the same expression, when written in terms 

^Of course, this is true provided that one imposes consistent renormahzation conditions. 
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Figure 6: Diagrams relevant for B^"^ and (1?+)^/^ in the absence of electro- 
magnetic corrections and isospin breaking effects. 

of functional integrals over quark and gluon fields. In the absence of 
electromagnetic corrections and isospin breaking effects, the same two 
diagrams as in fig. |^ contribute to both the decay amplitudes. In the 
presence of 5'f/(3)-breaking effects and of the electromagnetic interac- 
tion, which has to be included for consistency with the calculation of 
the Wilson coefficients, Bg^'^ is, however, different from (i?"*")^/^. We 
are not aware of any non-perturbative estimate of the difference of the 
S-parameters in the presence of the electromagnetic corrections. One 
may argue that the difference is small, and we have ignored it in our 
analysis. 

3. In ref. [|ll[], we have used the same value of Bg^^ as in the present 
analysis, Bg^'^{fi = 2 GeV) = 0.62, see table |^. This value should be 
compared to Bg^'^{fi = 2 GeV) ~ 0.44, which was used in ref. |TH]. Our 



estimate of Bg^'^ is derived from the theoretical value of Bk] in ref. [^], 
they evaluated Bg from the experimental value of Re A2 instead, by 
assuming SU{3) symmetry and by neglecting electromagnetic effects. 
It is well known that, from the AI = 3/2 amplitude, and using SU{3) 
flavour symmetry, one obtains a rather low value of the S^-parameter, 
and hence of Bg = Bx- This seems to us the origin of the difference. 
Given the uncertainties in the extraction of the S-parameters, and the 
electromagnetic and S'f/(3)-breaking effects discussed above, we think 
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that the two different values of -Bg^^ are both theoretically acceptable^. 

4. In the absence of a lattice calculation, the values of assumed in 
all our analyses have been guessed taking into account that the cor- 
responding matrix elements vanish in the vacuum insertion approxi- 
mation. Following ref. we have normalized the matrix elements 



of Ql, to those of Q^il^hj writing {QlM) = (Qifif^)) ■ 

The range of values used by us, given in table ^, has been criticized 
in ref. |]T3[. Their argument is the following. By matching at thresh- 
old {fi = rric) the effective Hamiltonian, valid for a scale /i > rric, 
to the effective Hamiltonian, valid for fi < rric, and by assuming the 
positivity of (O^.), they demonstrated that -Bi2(/^ = ^c) = 0. The 
values of -Bi2(/^ > ^c) were then evaluated by using the renormaliza- 
tion group evolution of the corresponding operators; they were found 
to be much smaller than the upper values used in our case. We be- 
lieve that this demonstration is not valid because an important point 
was missed. The effective Hamiltonian, valid for /i < rric, cannot be 
used to describe the physics at threshold. This is due to the presence 
of operators of dimension higher than 6, which are usually neglected 
at /i <^ rric, because their contribution decreases as {n/rric)^, with k a 
positive integer^. Indeed, for n <C rric, the effect of these higher dimen- 
sional operators enters only in the matching conditions; at threshold, 
however, all the higher dimensional operators become equally impor- 
tant and their contribution should be included. This turns out to be 
impossible in practice. For this reason, we believe that our guess for 
BI21 although with no theoretical foundation, is an acceptable one. 



8 Theoretical uncertainties 

In this section we briefly summarize the main sources of uncertainty present 
in the determination of e'/e. Since many aspects of the calculation have 

^^In ref. we simply wrote Bk — . This led the authors of ref. to believe 
that we had used Bl'^{^i = 2 GeV) - 0.8. 

^■^Indeed one may wonder whether a region of /i where it is possible to neglect terms of 
0(/i/mc) and at the same time to compute the Wilson coefficients in perturbation theory 
really does exist. 
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already been discussed in our previous papers on the same subject [TT 
(see also ref. |jl5[) we only add here a few considerations that may be useful 
to the reader. 



Theoretical predictions of weak decay amplitudes rely on the perturbative 
calculation of the Wilson coefficients and on the non-perturbative evaluation 
of the corresponding matrix elements. Even though the two aspects are 
intimately connected, we discuss them separately. 



8.1 The Wilson coefficients 

The main uncertainties in the evaluation of the Wilson coefficients are due to 
higher order corrections, which will probably remain unknown for a long time, 
to the value of Aqcd, which is taken from experiments, and to the choice of 
the expansion parameter (i.e. a^) and of the renormalization scheme. The 
latter two are also related to the truncation of the perturbative expansion 
and can be reduced by increasing the scale n at which the operators are 
renormalized. It should be noticed that the results found in HV and NDR, 
for a top quark mass of ~ 175 GeV, differ by about 22%. In principle, 
the differences of the Wilson coefficients in the two different renormalization 
schemes are compensated, at the NLO, by the corresponding differences in 
the renormalized operators, and hence in their matrix elements. The ma- 
trix Ar connecting the two different schemes (cf. sec. ^ and refs. 
has, however, quite large elements. This, combined with the cancellations 
occurring with a heavy top mass, is at the origin of the sizeable difference 
(although of 0{al)) found between the values of e'/e in HV and NDR. 



The larger differences observed in ref. [15| between HV and NDR are 
due to the fact that the authors changed the coefficient functions without 
modifying accordingly all the corresponing 5-parameters, in particular they 
kept the values of Bq and -Bg^^ fixed. 



8.2 Matrix elements 

The main limitation to an accurate theoretical prediction of e'/e does not 
come, however, from the calculation of the Wilson coefficients, but it is due 
to the evaluation of the operator matrix elements. We have chosen to work 
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in the framework of lattice QCD for several reasons. Firstly, the matching 
between the full theory and the effective Hamiltonian can consistently be 
done at the NLO (or at any higher order). This is not the case of other 
methods, like for example the l/A^^c-expansion. Secondly, the calculation 
of the hadronic matrix elements of the relevant operators is based on the 
fundamental theory without extra-assumptions: there are no free parameters 
besides the quark masses and the value of the lattice spacing, both of which 
are fixed by hadron spectroscopy. Thirdly, the typical scale at which the 
operators are renormalized is of the order of the inverse lattice spacing (see 
sec. which, in current lattice simulations, is relatively high, 1/a ~ 2-4 
GeV. A large renormalization scale, /x ~ 1/a, reduces the systematic error 
due to the presence of higher-order terms in the matching procedure, as can 
be seen from the fact that the NLO corrections are smaller at larger values of 
fi []TT], |T5|. Finally, the systematic errors present in lattice calculations (finite 
volume effects, discretization errors, "quenching", etc.) can be studied and 
reduced in time, according to the progress realized with the advent of more 
powerful computers. In spite of the many advantages offered by the lattice 
approach, the errors in the theoretical predictions are still rather large. This 
is due to a variety of reasons: 



Almost all the calculations have so far been done in the so-called 
"quenched" approximation, which consists in neglecting the effects of 
the virtual quark loops in the numerical simulation. This approxi- 
mation is a necessity imposed by the present limitation in computer 
power. The effects of the quenched approximation on the different 
physical quantities is not known a priori. In ref. |3l|, it has been 



shown that the value of the kaon 5-parameter Bk, obtained in a full, 
"unquenched" QCD calculation with staggered fermions, does not dif- 
fer appreciably from the corresponding one obtained in the quenched 
case. It should be noticed, however, that the mass of the quarks at 
which the calculation of ref. ^T[| was performed is rather large, thus 
reducing the effects of the quark loops (the quenched approximation 
can be obtained by sending to infinity the mass of the quarks in the 
loops) . 

Renormalized lattice operators must correspond, in the limit of infi- 
nite lattice cut-off, to finite chiral covariant operators, which obey the 



38 



same renormalization conditions as in the continuum. With few ex- 
ceptions, lattice perturbation theory is usually used to evaluate the 
renormalization constants of lattice operators. The problem of mixing 
with lower dimensional operators requires, however, a non-perturbative 
subtraction of all power divergences [^-[^. Apart from this 



special case, the use of perturbation theory, in the computation of 
multiplicative constants and adimensional mixing coefficients, is well 
justified, provided that the lattice spacing a is sufficiently small, i.e. 

^ Aqcd- However, there is evidence of failures of lattice per- 
turbation theory, resulting from the use of the bare lattice coupling 
constant a'"** = Qq/Att, as an expansion parameter. This failure has 
been attributed to the presence of "tadpole" diagrams, which appear 
only in the lattice regularization and are likely to give large corrections 
in higher order perturbation theory. Several solutions to this problem 
have been proposed so far [^, |55|, which could hopefully reduce the 
systematic error induced by the poor convergence of lattice perturba- 
tion theory. It remains to be shown that these methods work in the 
case of the four-fermion operators relevant for e'/ e. 

In most of the lattice simulations, the {TT\Oi\K) matrix elements, in- 
stead of the more appropriate {7ni\0i\K) ones, are indeed computed. 
The K ^ 71 amplitudes are then related to the physical K ^ mi ones 
using current algebra, at lowest order in the chiral expansion. This 
certainly introduces an error, which on a single matrix element is prob- 
ably less than 30%. Also in this for the isospin breaking effect 
mentioned before, the systematic error is amplified by the large cancel- 
lations occurring in the final result. 

The crucial 5-parameter Bq has been computed so far (together with 
B^) only by one group using "staggered" lattice fermions. 



All the attempts to compute -85^6 using different lattice formulations of 
QCD, typically using "Wilson fermions" , have failed so far. Morevover, 
the same techniques applied to a similar problem, the A J = 1/2 
K —> 7171 amphtude, were completely unsuccessful. It would be more 
reassuring to have further confirmations of the values of i?-parameters 



computed in refs. [23, 24 
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9 Conclusions 



In this paper we have discussed several theoretical issues, which are relevant 
to the prediction of e'/e and many other weak decays. To our knowledge, the 
discussion of the matching of the lattice operators to the continuum effective 
Hamiltonian is new. Sections ^ and ^ contain several useful formulae that can 
be used in many different applications. With the operator matrix elements 
from the lattice, one can compute the weak decays amplitudes, using the 
Wilson coefficients taken from table 4 (5) or 6 (7), with almost no extra 
effort. The RI renormalization scheme discussed in this paper will become 
very useful if the non-perturbative method proposed in ref. will succeed 
in the case of four-fermion operators. 

We have also presented an upgraded theoretical analysis of e'/e in two 
different renormalization schemes, NDR and HV. Let us summarize the 
main results of the present study. Using the information coming from e and 
Xd, and taking nit = (174 ± 17) GeV and fsB^^^ = (200 ± 40) MeV, we 
obtain constraints on the angles of the unitary triangle 

cos 5 = 0.47 ±0.32, 

sin 2/3 = 0.65 ±0.12. (70) 
Moreover, for cos 5 > 0, we predict 

= 0.75 ±0.23. (71) 

Finally, we have obtained 

e'/e = (3.1 ±2.5) X 10"^ (72) 

as our best estimate. Isospin-breaking effects have been found quite im- 
portant, because of the large value of m^, and we believe that they deserve 
further theoretical investigation. 
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